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Abstract 

We consider three optomechanical cavities coupled irreversibly as well as reversibly to each other 
and explore entanglement between the different optical and mechanical modes in the collective 
system. Each cavity is driven by a coherent field and the optomechanical coupling in each cavity is 
treated in the linearised regime. The composite system exhibits robust intercavity photon-phonon 
entanglement, in bipartite form, well into the steady state. 
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I. INTRODUCTION 



Quantum cavity optomechanics capitalizes on the radiation pressure of light exerted on 
a mechanical degree of freedom incorporated in confined high finesse cavities [IH3] to en- 
able coupling between an optical and mechanical resonator. This currently active area of 
quantum physics is witnessing accelerated theoretical development p-[TO] as well as ex- 
perimental achievements jlTHTS] . Typically the cavity is driven to a steady state and the 
optomechanical coupling is linearised around the steady state amplitude. More recently, vac- 
uum optomechanical coupling was studied (TB] . Hence optomechanics serves as an excellent 
test bed for fundamental experiments at the quantum-classical boundary, paving innovative 
ways towards controlling the mutual interaction between light and the mechanical motion of 
mesoscopic objects. Sideband cooling to the absolute ground state of mechanical resonator 
systems [T7H19] turned out to be quite a challenge experimentally, due to the detrimental but 
unavoidable thermal coupling of the resonators to their environments. Nevertheless these 
challenges have now been overcome, with different experimental groups recently having suc- 
cessfully achieved the ground state of mechanical oscillators in cavity electromechanical and 
optomechanical systems [20H23] . Hence we are now in the era of ground state mechanical 
systems and finally, the long anticipated quantum mechanical control of macroscopic objects 
is in sight. See [2^ for an overview and outlook. 

Akram et al. [25j recently showed how a discrete exchange of a non classical state, the 
single photon, can take place between the optical and mechanical modes of an optomechan- 
ical system, providing a feasible route towards quantum state transfer. In this paper we 
are interested in the dynamics of coupled optomechanical arrays. The mechanical and op- 
tical modes within one optomechanical system have been shown to exhibit a considerable 
degree of entanglement [261429] . The entanglement between the output optical fields of a 
trapped-mirror-system has also been shown in [SOj . A remarkable feature of optomechanical 
entanglement is that it can exist at temperatures above the absolute ground state. In pi] . 
Mazzola and Paternostro consider the entanglement between a pair of optomechanical cavi- 
ties, in the linearised regime, that arise when each cavity is driven by one of the twin beams 
generated by a source of spontaneous parametric down conversion. The optomechanical 
systems become entangled due to quantum correlations in the light sources driving each 
cavity. Here we explore the entanglement between optical and mechanical modes of three 
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different optomechanical cavities where they are coupled via their optical ports in either a 
reversible or irreversible (cascaded) manner. Another approach to coupled optomechanical 
arrays would be to couple the mechanical resonators rather than the optical resonators • 
A system of two coupled microwave cavities each containing a mechanical element has also 
been considered by Heinrich and Marquardt [33]. 

II. THEORY AND MODEL 

The basic idea in quantum cavity optomechanics is to induce a reversible coupling be- 
tween an optical and mechanical resonator. Typically the interaction arises from the ra- 
diation pressure of light. The usual set up is modelled around an optical cavity whose 
resonance frequency is altered by the displacement of some mechanical resonator modelled 
as a harmonic oscillator. A shift in the resonance frequency of the optical cavity changes 
the circulating power and thus changes the radiation pressure on the mechanical resonator, 
yielding the optomechanical coupling which gives rise to a plethora of effects depending on 
how the various parameters and configurations are manipulated in the system. For example, 
the cavity can be driven to a steady state amplitude and the vacuum interaction linearised 
around this amplitude. This gives a coupling that is quadratic in the amplitude of the 
optical and mechanical resonator. The goal in recent experiments has been to push the 
macroscopic mechanical elements of optomechanical systems towards the quantum limit by 
various passive and active cooling protocols. 

In this paper we consider a composite system of a chain of three identical cavities where 
the optical and mechanical resonators are quadratically coupled to each other. This can be 
accomplished in two different ways. Each cavity (1,2,3) is an optical Fabry-Perot cavity in 
which one of the mirrors is subject to a harmonic restoring force and can thus move due to 
radiation pressure, or with both mirrors fixed and an optically levitating nanosphere trapped 
within the optical cavity. Recently, optically levitating nanospheres have been proposed as 
new optomechanical systems with much higher Q factors due to their natural isolation from 
a thermal contact point [OHn] which is the most detrimental obstacle in achieving the pure 
ground state of the mechanical oscillator via optomechanical cooling. It is expected that 
these particular optomechanical systems will be much more efficiently cooled to their ground 
states as even at room temperature they are excellently well isolated [6]. 
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The mechanical resonator in either model has a frequency Wmi,2,3 damping rate ^1,2,3 
while the optical cavity has a resonance frequency ^ 3 and is strongly driven with a co- 
herent pump field at frequency u^-^ ^ 3 . Further, the trio of optomechanical cavities can 
be coupled together in two different configurations: via an irreversible coupling or a re- 
versible coupling configuration as illustrated in Fig. [T} 



circulator circulator 




Reversible coupling 



FIG. 1: Model: The three optomechanical cavities can be coupled either irreversibly via forward 
feed tunneling using circulators or reversibly via nearest neighbours. 

The irreversible coupling is a forward feed method modelled as the cascaded systems 
approach |3l]. Reversible coupling requires that the cavities be evanescently coupled and 
thus must be close to each other, however the irreversible (cascaded) coupling does not 
require this, but wiU require circulators be placed between the optical cavities to ensure 
irreversibility. 

The Hamiltonian for the composite system is given as 

3 

H = h'^Aattltta + UJmjiba + Qaitta + + (l) 

a=l 

where a/b and a^fo^ are the annihilation and creation operators respectively of each opti- 
cal/mechanical mode and g is the optomechanical coupling rate within each cavity. Here 
we have considered the shifted reference frame where the amplitude of the cavity has been 
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displaced in the steady state due to the driving laser. Hence g is the effective couphng 
strength proportional to the steady state amplitude of the cavity field, due to linearisation 
of the radiation pressure force. A„ = uj^^ — ujl^, is the detuning with respect to the driving 
laser field with a = 1,2,3. The master equation for the composite system in the interaction 
picture with respect to the driving laser 2,3 given as, 

a=l 

revH 

+ 

where the reversible and irreversible dynamics are given respectively by 

J^rev = -«Xl2[ai4 + «1«2, p] " «X23[a24 + 4«2, p] (3) 

Cirr = y/niH2{[aip, a^] + [a2, pa\]) + y/KiK3{[aip, 4] + [as, pa\]) (4) 
+^/K^{[a2P, 4] + [«3, pal]) 

We assume that the coupling is either reversible or irreversible but not both together. Here 
Kj is the line width of each cavity-j, while Xkj is an arbitrary coupling strength between the 
cavities k and j for the reversible coupling configuration and the damping superoperator V 
is defined by 

V[A]p = ApA^ - liA^Ap + pA^A) (5) 

As the cavities are driven by an external field, the coupling between the cavities is induced 
after each system attains steady state. 

III. COUPLED OPTOMECHANICAL ARRAYS 
A. ForwEird feed coupling 

We first consider the forward feed coupling configuration between the three optome- 
chanical cavities. Here the cavities are coupled via a unidirectional coupling only. In our 
calculations, we have neglected time delay between the cavities. Effectively, we can consider 
cavity 1 as the " main source" cavity while cavity 3 is only a receiver cavity. Cavity 2 on the 
other hand receives photons from cavity 1 as well as drives cavity 3. The topology of our set 
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up in this coupUng configuration is further illustrated in Fig. [2j Hence we have reversible 



interactions between the optomechanical branches of the chain but irreversible interactions 
between the optical ports coupling the individual optomechanical units. This is quite a 
different configuration to previous work on coupled oscillator arrays [35j . 



FIG. 2: Topology of chain of three optomechanical systems for forward feed (irreversible) coupling 

As has been shown previously ^29j, intracavity photon-phonon entanglement is present 
within each optomechanical unit. However here we are interested in the presence of entan- 
glement between the optical and mechanical resonators in distinct optomechanical cavities. 
Our composite system under study is in the Gaussian state as it starts from the vacuum, 
and the equations of motion are linear, we will have all Gaussian states in the system. 
Consequently, to quantify the entanglement, we use the logarithmic negativity measure for 
Gaussian states formulated by Vidal [36] . The log negativity between two states is expressed 
in terms of the entries of their covariance matrix, 7 which is a 4 x 4 matrix given as 



The matrices 7a,b arise from position and momentum quadratures of the optical and 
mechanical modes respectively, while 7c is a result of cross terms between the optical and 
mechanical position and momentum quadratures. The log negativity is then calculated as 





(6) 





otherwise. 



(7) 
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FIG. 3: Forward feed coupling 

The log negativity between the intercavity phonons of the composite system of the three optome- 
chanical units (61,62,^3) vs. time for ki_2,3 = 0.1, ^1,2,3 = 0.001, Ai_2, 3 = ^1,2,3 = 40,51,2,3 = 
0.5, n = 0. 

where the function, / : C(4) — t- SR"*" is defined as 



/(7) = rA,B,c - ^JTXb,c - l7l (8) 

such that 

p I I + I 7b I I I /ON 

^A,B,C = ^ I 7c I (9) 

Hence in order to determine the matrix elements of each of the 2x2 matrix 7yi,B,c we 
calculate the second order moments from the master equation of the composite system. 
Inserting the relevant second order moments into Eq. [7] we are now able to calculate the 
entanglement between any two modes of the system. 

Fig. [3] shows the bipartite entanglement between all the mechanical mirrors &i,&2 and 63 
of the composite system for the feed forward coupling between the optical ports. We note 
the intercavity phonon-phonon entanglement exhibited by this system lasts well into the 
steady state. 

Fig. |4] shows the entanglement between all possible pairs of the optical and mechanical 
modes of different cavities, i.e. intercavity photon-phonon entanglement. Clearly all the 
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FIG. 4: Forward feed coupling 
The log negativity between all possible selections of intercavity photons and phonons of the com- 
posite system of three optomechanical units vs. time for: ki^2 = 0.1, ^1^2 = 0.001, Ai_2,3 = wi,2,3 = 
40,(71,2,3 = 0.5, n = 0. 

six modes of the system are shown to exhibit two mode bipartite entanglement within 
themselves. We note the symmetry for pairs of oscillators ai,bj and aj,bi. This is due to 
the symmetry in the Hamiltonian, Eq. [Tjfor cavities with identical parameters cu^i and gi. 

We have found the entanglement to be maximum when the detuning is resonant with 
the mechanical frequency, i.e. A = Um, hence the system is on the red sideband. We would 
like to remind the reader, that the full interaction Hamiltonian has been considered here, so 
that the squeezing terms dibi and dj^&J also contribute towards the dynamics of the system 
and in fact are responsible for the presence of the entanglement within each optomechanical 
unit. The traditional assumption, however would be to achieve strong entanglement over 
the blue sideband when A = —Um- But the amount of entanglement that could possibly be 
generated within each optomechanical unit is limited by the stability condition g < y/KjI/2 
arising in the Rotating Wave Approximation (RWA) limit, A = —Um » k,, [37]. However 
for optomechanical systems the mechanical damping rate /i is much less than the optical 
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damping rate, k, consequently limiting the coupling strength g within each optomechanical 
unit to very small values {g^O.Ol). As a result, for such weak coupling parameters, the 
intracavity entanglement does not transfer to intercavity modes on the blue sideband. 

On the red sideband however, while the beam splitter interaction dominates over the 
squeezing interaction, stability conditions for A = Um, allow evolution of the system under 
the strong coupling regime, g < uj'^ + , taking the system away from the RWA 
limit [33. This results in an exchange of excitations on the optomechanical branches of the 
chain of oscillators, consequently distributing the intracavity entanglement between oscil- 
lators of different cavities due to the forward feed coupling between the optical ports. It 
is interesting to see that photons and phonons of different cavities are robustly entangled 
without having been in direct contact with each other. Hence the forward feed coupling 
between the optomechanical units facilitates the distribution of intracavity entanglement 
over the intercavity modes. 

Applying the relevant inequalities from Ref. [3H], we determine that the multimode dis- 
tribution of the entanglement we are observing in these results does not exhibit any genuine 
tripartite distribution. Hence we conclude that here our composite system of a chain of 
three optomechanical systems forms a class of entangled states of the optical and mechani- 
cal modes. 

B. Reversible coupling 

We now investigate the presence and possible distribution of entanglement between inter 
cavity modes in the reversible coupling configuration. Here the three cavities are evanes- 
cently coupled with nearest neighbours with an arbitrary coupling strength Xij given in 
Eq. |3} The topology we are considering here thus allows for a reversible exchange of optical 
excitations between the cavities as shown in Fig. |5| 

In Fig. [6] we have the entanglement between the intercavity mechanical modes versus 
time. Clearly here also the mirrors of different optomechanical cavities are entangled for a 
long duration in the absence of a direct coupling induced between them. 

Fig. [7] shows steady state intercavity photon-phonon entanglement between all two mode 
bipartite selections of the composite system of coupled optomechanical units. Again the 
system is limited by different stability conditions on the blue and red sidebands as discussed 
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FIG. 5: Topology of chain of three optomechanical systems for reversible coupling 
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FIG. 6: Reversible coupling 

The log negativity between the inter cavity phonons of the composite system of the three op- 
tomechanical units (61,62,^3) vs. time for ki,2,3 = 0.1,/xi,2,3 = 0.001, Ai^2,3 = ^1,2,3 = 20, 31,2,3 = 
1,X12 = X23 = 0.5, n = 0. 



in section III A so that we find the entanglement is distributed over intercavity modes only 
for the red sideband. Hence the intercavity entanglement is found to exist only on the red 
sideband of the system, as the beam splitter part of the Hamiltonian induces the distribution 
of the intracavity optomechanical entanglement over the cross cavity modes in the strong 
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FIG. 7: Reversible coupling 

The log negativity between all possible selections of inter cavity photons and phonons of the 
composite system of three optomechanical units vs. time for: ^1,2, 3 = 0.1,/ii^2,3 = 0.001, Ai^2, 3 = 
wi,2,3 = 20,51,2,3 = 1,X12 = X23 = 0.5, n = 0. 

coupling regime, satisfying the stability conditions of the system. 
IV. SUMMARY 

We have found that considerable robust steady state intercavity pair-wise entanglement 
exists between all the six modes (01,2,3, &i, 2.3) of a chain of three optomechanical systems 
when they are coupled both irreversibly via a feed forward as well as reversibly, via a 
nearest neighbour evanscent coupling. One would expect that the entanglement between 
the three cavity modes would exhibit some form of tripartite entanglement however, we find 
that applying the relevant inequalities from Ref . [38] , the system of coupled optomechanical 
arrays does not exhibit genuine tripartite entanglement. Hence the system evolves such 
that pairs of oscillators are always accessible in two mode bipartite selections. Further, 
the topology of the chain of oscillators considered in this paper comprising of reversible 
optomechanical coupling between light and mechanics, alongside irreversible forward feed 
tunneling between optical ports, is quite distinct from earlier results on coupled oscillators. 
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We note that the oscillators are entangled for a longer duration in the irreversible coupling 
configuration as compared to the reversible case. The reversible case, based on nearest 
neighbour evanescent coupling, will require the cavities to be in close proximity, however 
the feedforward case is not so restricted. In this paper we have neglected delay effects in 
the feedforward case which means that delays must be less than the cavity decay times, 
which sets the time scale for the dynamics of the optomechanical array. Even with this 
assumption, the feedforward case enables the optical cavities to be much further apart — 
up to one metre — than the reversible coupling case. Coupled optomechanical arrays could 
have wide applications as components of quantum repeaters and quantum memories required 
in quantum information processing. 
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